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Final  Report 

Millimeter- Wave  Acoustic  Transducers 
June  1,  1987  -  September  30,  1989 


I.  Introduction 

This  project  was  focused  on  a  program  to  realize  efficient  systems  for  generating  and 
detecting  coherent  phonons  with  millimeter  wavelengths.  The  primary  task  was  the  realization  of 
transducers  that  are  capable  of  converting  electromagnetic  energy  to  acoustic  energy  in  this 
frequency  region. 

In  order  to  increase  the  operating  frequency,  the  components  of  the  acoustic  transducer  must 
be  redesigned.  Instead  of  the  single  layer  ZnO  acoustic  transducers  used  below  10  GHz, 
multilayer  ZnO  acoustic  transducers  must  be  used  to  generate  sufficient  power  at  frequencies  above 
24  GHz.  These  multilayer  ZnO  transducers  have  been  used  to  generate  approximately  28  GHz 
sound  with  two-way  conversion  losses  as  low  as  27  dB.  Multilayer  transducers  have  also  been 
fabricated  for  operation  at  frequencies  as  high  as  100  GHz,  where  the  two-way  conversion  loss  is 
approximately  50  dB.  These  ZnO  multilayer  transducers  represent  the  highest  operating  frequency 
of  any  ZnO  transducer  and  the  highest  conversion  efficiencies  for  any  type  of  transducer  in  this 
frequency  range.  The  theoretical  aspects  of  multilayer  transducers  were  explored  along  with 
experimental  verification  of  some  theoretical  predictions.  The  possibility  of  generating  sound  at 
frequencies  as  high  as  300  GHz  has  been  considered. 

A  piezoelectric  mechanical  scanning  system  was  designed  for  operation  at  temperatures  as 
low  as  100  milliKelvin.  The  scanner  is  able  to  overcome  the  loss  in  sensitivity  of  the  piezoelectric 
elements  at  low  temperatures  by  using  lever  arms  to  magnify  the  piezoelectric  motion.  The 
resolution,  frequency  characteristics,  and  heat  dissipation  properties  of  the  scanner  are  presented. 

In  association  with  all  of  this  we  have  continued  to  investigate  the  nature  of  acoustic  wave 
propagation  in  liquid  helium.  This  liquid  is  suitable  for  propagating  sound  waves  at  these  ultrahigh 
frequencies  and  our  goal  was  to  study  the  Unear  and  nonhnear  properties  of  this  medium. 
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As  a  final  effort  near  the  end  of  this  period  we  co-operated  with  research  workers  at  the 
Aerospace  Corporation  in  El  Segundo  to  construct  a  scanning  acoustic  microscope  operating  at 
15.3  GHz  with  a  resolution  of  150  A. 

II.  Millimeter-wave  Acoustic  Transducers 

Previous  generations  of  water  and  cryogenic  acoustic  microscopes  used  thin  film 
piezoelectric  transducers  to  generate  and  detect  the  acoustic  signal  used  for  imaging.  These 
transducers  usually  consist  of  a  vacuum  deposited  thin  film  of  piezoelectric  material  sandwiched 
between  two  metal  electrodes.  An  rf-voltage  is  applied  across  the  electrodes  and  acoustic  waves 
are  generated  in  the  piezoelectric  material.  In  a  reciprocal  manner,  acoustic  waves  entering  the 
piezoelectric  material  generate  a  voltage  across  the  electrodes  which  can  be  sensed  and  recorded. 
Such  transducers  generally  operate  at  thickness-mode  resonance  and  are  probably  the  most  efficient 
generators  and  detectors  of  acoustic  waves  for  frequencies  ranging  from  a  few  hundred  Megahertz 
up  to  10  GHz.1’2  Thin  film  transducers  are  capable  of  generating  and  detecting  both  shear  and 
longitudinal  waves  depending  on  the  orientation  of  the  piezoelectric  material  with  respect  to  the 
applied  electric  field.  A  large  variety  of  materials  have  been  used  to  generate  microwave  frequency 
acoustic  waves  including  ZnO,  CdS,  and  LiNb03.  In  transducers  reported  here,  ZnO  was  chosen 
because  of  its  high  coupling  coefficients  and  ease  of  vacuum  deposition.  At  frequencies  above 
10  GHz,  the  performance  of  these  simple  acoustic  transducers  is  impaired  by  two  factors.  The 
conversion  efficiency  typically  becomes  quite  small  for  reasonably  sized  transducers  and  the  failure 
rate  from  fabrication  defects  increases  with  the  decreasing  thickness  of  the  transducer.  In  order  to 
extend  the  operating  frequency  of  the  acoustic  microscope,  a  new  approach  to  the  transducer 
construction  was  needed.  This  chapter  describes  the  design,  fabrication  and  evaluation  of 
multilayer  ZnO  acoustic  transducers  in  the  frequency  range  of  25-100  GHz  and  explores  the 
possibility  of  using  multilayer  ZnO  transducers  to  generate  sound  at  frequencies  above  100  GHz. 

1.  Rational  for  Multilayer  Transducers 

The  conversion  loss  of  a  thickness-mode  resonant  transducer  can  be  separated  into  two  parts: 
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CL(dB)=CLac(dB)+  CLel(dB) 


(2.1) 


where  CLac  is  the  part  due  to  the  acoustic  losses  in  the  system  and  CLei  is  the  part  due  to  the 
electrical  losses  in  the  system.  Eq.  (2.1)  can  be  simplified  for  many  microwave  transducers.  In 
most  microwave  applications,  the  thicknesses  of  the  electrodes  are  small  compared  to  an  acoustic 
wavelength  and  their  effect  on  the  conversion  loss  can  be  neglected.  If  there  are  no  acoustic 
matching  layers  in  the  system,  the  electrical  losses  dominate  the  acoustic  losses  since  most 
thin-film  materials  have  high  acoustic  Q's.  The  electrical  loss  can  be  characterized  by:1 
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where  Rs  is  the  source  impedance,  Ra  is  the  real  part  of  the  transducer  impedance,  Xa  is  the 
imaginary  part  of  the  transducer  impedance.  R,*  and  X^  represent  the  real  part  and  imaginary  part 
of  the  stray  series  impedance  that  usually  arises  from  the  physical  connection  of  the  transducer  to 
the  source.  For  most  applications,  Rs  is  typically  50Q,  Rse  and  Ra  are  on  the  order  of  IQ,  and  Xse 
is  usually  small  until  very  high  frequencies.  If  Xa  «  Rs,  Eq.  (2.2)  can  be  approximated  by: 

CLel(dB)=  lOlog 

Because  Rs  is  set  by  the  equipment,  the  major  factor  affecting  the  conversion  loss  is  Ra. 

If  the  acoustic  losses  are  small  and  the  thickness  of  the  metal  electrodes  is  much  less  than 
the  sound  wavelength,  the  real  part  of  the  transducer  impedance  at  resonance  can  be  approximated 
as: 
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where  Jq  is  the  piezoelectric  coupling  constant  of  the  material,  r^  is  the  ratio  of  the  mechanical 
impedances  of  the  thin  film  and  the  substrate,  (Oo  is  the  angular  frequency  at  the  half-wave 
resonance  and,  Co  is  the  clamped  capacitance  of  the  transducer.  The  clamped  capacitance  of  the 
transducer  is  inversely  proportional  to  the  half-wave  thickness  and  proportional  to  the  resonant 
frequency.  Thus,  for  a  fixed  r^,  the  radiation  resistance  decreases  as  the  resonant  frequency 
squared.  The  result  of  Eqs.  (2.3)  and  (2.4)  is  that  if  one  simply  tries  to  scale  the  center  frequency 
of  the  transducer  from  10  GHz  to  100  GHz,  the  two-way  untuned  conversion  loss  increases  by  40 
dB  as  can  be  seen  in  Fig.  2.1. 

Several  methods  have  been  proposed  to  increase  the  conversion  efficiency  at  high  gigahertz 
frequencies.  The  conversion  loss  can  be  improved  at  the  expense  of  reduced  bandwidth  by 
electrically  tuning  the  transducer  to  cancel  the  impedance  mismatch  between  the  source  and 
transducer.  If  the  impedance  mismatch  were  able  to  be  perfecdy  compensated,  a  difficult  task  at 
millimeter- wave  frequencies,  the  tuned  conversion  loss  would  be: 


CLel(dB)=101og 
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At  millimeter- wave  frequencies,  two  factors  dominate  the  series  resistance.  Most  methods  for 
attaching  the  transducer  to  the  rf  source  employ  thermo-compression  bonding  to  the  top  electrode 
which  generates  some  resistance  at  the  bond  interface.  The  other  major  factor  affecting  the  series 
resistance  is  the  small  skin  depth  of  the  rf  signal  at  high  frequencies  which  leads  to  high  surface 
resistivities.  These  two  factors  usually  place  a  lower  limit  on  the  series  resistance  of  about  1  £T 
The  impedance  of  a  single  layer  transducer  can  be  estimated  from  Eq.  (2.4).  A  60  |im  radius  ZnO 
transducer  operating  at  a  center  frequency  of  100  GHz  on  a  sapphire  substrate  would  have  an 
impedance  of  less  than  0.01£I  Even  in  such  a  best  case  example,  the  tuned  one-way  conversion 
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Conversion  Loss  (dB) 


Figure  2.1  Two-way  conversion  loss  of  a  thin-film  transducer  as  a  function  of  frequency. 
Transducers  radius  =  60  |im,  sapphire  substrate. 
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loss  would  still  be  approximately  20  dB  making  tuning  alone  an  impractical  approach  to  the 
problem. 

Another  possible  solution  is  to  increase  the  transducer's  impedance  by  decreasing  the 
transducer  radius.  In  the  high  frequency  water  microscope,  transducers  with  radii  as  small  as 
30  pm  have  been  used.3  In  the  helium  acoustic  microscope,  reducing  the  transducer  radius  would 
require  decreasing  the  lens  radius  to  maintain  uniform  aperture  illumination.3’4  In  addition  to  the 
difficulties  in  the  precise  machining  of  a  smaller  lens  opening,  the  overall  length  of  the  lens  would 
have  to  be  reduced  to  ensure  proper  illumination.3  To  overcome  the  small  reduced  path  length  in 
the  sapphire,  shorter  pulses  would  have  to  be  used  at  larger  bandwidths-leading  to  higher  noise 
levels  and  reduced  dynamic  range.  Because  of  the  numerous  detrimental  side  effects,  reducing  the 
transducer's  radius  as  the  sole  method  of  improving  the  conversion  efficiency  is  undesirable. 

The  second  problem  encountered  in  using  single  layer  acoustic  transducers  at  millimeter  wave 
frequencies  is  that  the  small  thickness  of  the  transducer  would  often  lead  to  shorting  between  the 
two  metal  electrodes  rendering  the  transducer  inoperative.  Such  shorting  was  observed  in 
transducers  fabricated  for  operation  at  8  GHz  where  the  ZnO  is  approximately  2000  A  thick.3 
Near  100  GHz,  the  typical  thickness  for  a  ZnO  transducer  would  be  a  few  hundred  Angstroms 
increasing  the  likelihood  of  pinholes.  The  occurrence  of  pin-hole  defects  can  be  reduced  by 
improving  the  sputtering  uniformity;  however,  at  millimeter  wave  frequencies  shorting  between  the 
electrodes  could  still  be  a  serious  problem. 

2.  Multilayer  Theory 

In  an  effort  to  overcome  the  high  conversion  loss  associated  with  high  frequency  thin  film 
transducers,  Shaw5  and  de  Klerk6  proposed  that  multilayer  transducers  be  used.  In  a  multilayer 
transducer,  piezoelectrically  active  layers  alternate  with  passive  or  piezoelectrically  less  active 
layers  to  form  an  acoustic  resonator.  For  maximum  conversion  efficiency,  the  thickness  of  each 
layer  should  be  one-half  an  acoustic  wavelength  for  that  material.  Early  work  concentrated  on  CdS 
as  the  active  material  and  SiO6  as  the  passive  material.  In  the  early  1970's,  the  advent  of 
superlattices  made  it  theoretically  possible  for  coherent  phonon  generation  with  frequencies  into  the 
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infrared.7-8  Other  authors9-10  have  suggested  that  a  nonpiezoelectric  semiconductor  with  periodic 
doping  in  the  pattern  n  doped-/ ntrinsic-p  doped-/ntrinsic-n-/-p-/...  could  be  used  to  generate  high 
frequency  coherent  phonons  via  Coulomb  forces  between  the  accumulated  charges  at  the 
interfaces. 

Recently,  Hadimioglu  et  aln  reported  a  new  fabrication  technique  for  multilayer 
transducers.  This  technique  exploits  the  fact  that  ZnO,  and  all  hexagonal  crystals,  exhibit  a 
maximum  in  the  longitudinal  wave  coupling  when  the  c-axis  is  oriented  parallel  to  the  applied 
electric  field  and  a  minimum  at  approximately  40  degrees  from  the  electric  field.12  Other 
researchers13  have  shown  that  "tilted"  ZnO  can  be  grown  at  an  oblique  angle  to  the  c-axis  by  tilting 
the  substrate  and  positioning  it  near  the  edge  of  the  plasma  in  the  deposition  station.  By  altering 
the  c-axis,  and  thus  the  longitudinal  wave  coupling,  a  multilayer  transducer  can  be  assembled  one 
layer  at  a  time  as  shown  in  Fig.  2.2.  The  ZnO  is  deposited  by  dc-magnetron  sputtering  using  a 
zinc  target  in  a  80%-20%  oxygen-argon  atmosphere.  The  basic  ZnO  deposition  system  used  in 
depositing  single  layer  ZnO  transducers  has  been  modified  so  that  the  substrate  can  be  positioned  at 
three  different  locations  as  shown  in  Fig.  2.3(a).  In  position  I,  the  substrates  are  parallel  to  the 
ZnO  target  and  centered  over  the  target  When  deposited  in  this  configuration,  the  ZnO  grows  in 
the  "normal"  orientation  yielding  a  high  longitudinal  piezoelectric  coupling  coefficient.  In  position 
n,  the  substrates  are  rotated  40°  with  respect  to  the  ZnO  target  normal.  The  center  of  the  substrate 
holder  is  also  moved  a  distance  d  =  1.5  inches  from  the  location  at  position  I.  ZnO  deposited  in 
this  "tilted"  position  yields  a  smaller  longitudinal  piezoelectric  coupling  constant  ZnO  deposited  in 
position  III  has  similar  characteristics  to  that  deposited  at  position  n.  Both  positions  are  used  to 
ensure  that  any  thickness  variations  in  the  direction  parallel  to  the  substrate  holder  are  minimized. 
The  substrate  holder  is  moved  between  the  positions  in  the  pattern  I-II-I-in-I-n...  until  the  desired 
number  layers  have  been  deposited.  As  can  be  seen  in  Fig.  2.3(b),  the  substrate  holder  does  not 
have  circular  symmetry .  This  anisotropy  can  lead  to  transducers  with  different  performance 
characteristics  and  center  frequencies  depending  on  the  substrate  location. 
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Ka-BAND  TRANSDUCER  GEOMETRY 


Figure  2.2  Ka-Band  transducer  geometry. 

Schematic  showing  the  basic  elements  of  a  multilayer  transducer. 
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holes  for  substrates 
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Figure  2.3  Multilayer  transducer  deposition  schematic. 

(a)  In  position  I  the  ZnO  is  deposited  in  the  "normal”  normal  orientation,  in  positions  II 
and  HI  the  ZnO  is  deposited  in  the  "tilted"  orientation,  (b)  The  "head  on"  view  of  the 
sample  holder  showing  the  lack  of  circular  symmetry. 
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The  method  described  above  offers  several  advantages  over  other  methods  proposed  to  date. 
Because  only  ZnO  is  used  in  the  transducer  fabrication,  the  entire  deposition  can  be  done  in  a 
single  pump-down  in  a  dedicated  vacuum  station,  whereas  the  CdS  multilayer  transducers  used 
other  components  for  the  passive  layers.  Because  the  transducers  can  be  deposited  in  a  single 
pump-down  of  the  vacuum  system,  the  chances  of  contamination  are  minimized  and  the 
repeatability  of  the  results  is  enhanced.  Two  component  multilayer  transducers  also  present 
additional  problems  regarding  lattice  mismatch  between  the  layers8  and  stresses  developed  by 
differences  in  the  thermal  contraction  of  the  layers  between  room  temperature  and  typical  cryogenic 
acoustic  microscopic  applications. 

In  order  to  get  a  feel  for  the  advantages  of  a  multilayer  transducer,  one  can  imagine  building 
an  ideal  multilayer  transducer.  This  transducer  would  be  composed  of  N  layers  of  an  active 
piezoelectric  material  (A;;  i  =  1...N)  and  N  layers  of  a  passive  material  (Pj ;  i  =  1...N)  with  the 
same  acoustic  impedance  and  electrodes  so  thin  that  their  effect  is  negligible.  Each  of  the  layers 
would  be  one-half  an  acoustic  wavelength  thick  and  the  acoustic  attenuation  would  be  zero.  Such  a 
transducer  is  shown  in  Fig.  2.4(a).  Because  the  electromagnetic  radiation  pattern  is  of  nearly 
constant  phase  between  the  electrodes,  each  of  the  active  layers  will  have  the  same  response  as 
shown  in  Fig.  2.4(b).  The  acoustic  wave  entering  the  second  active  layer  will  have  the  same  phase 
as  the  acoustic  wave  generated  in  that  layer  and  the  amplitudes  will  add  as  scalars.  This  process 
continues  down  the  length  of  the  transducer  where  the  final  amplitude  will  be  N  times  the 
amplitude  generated  by  a  single  period.  Fig.  2.4(c).  Therefore,  the  acoustic  power  delivered  to  the 
substrate  is: 


P  «  A2  «  N 2 

a c  a c 


(2.6) 


where  Pgc  is  the  acoustic  power,  and  Aac  is  the  acoustic  wave  amplitude.  This  simple  analysis 
show  the  major  advantage  of  the  multilayer  transducer,  namely  that  the  conversion  efficiency  of  the 
transducer  is  proportional  to  the  number  of  layers  squared. 
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X 


Figure  2.4  Acoustic  wave  amplification  in  a  multilayer  transducer. 

(a)  A  transducer  containing  N  active  layers  (Aj)  and  N  passive  layers  (P,).  (b)  Propagation 
of  acoustic  wave  of  amplitude  Ao  generated  at  each  interface,  (c)  Amplification  resulting 
from  all  acoustic  waves  adding  in  phase. 
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A  more  complete  analysis  of  the  acoustic  power  generated  by  the  multilayer  transducer  can  be 
found  in  the  analysis  of  Yang  et  al  ?  This  analysis  yields: 

AgcfN.)  =  M(Fj  +  F2)  (2.7) 


where: 


M  = 


j  (k-q)Nt 

i  -e 


j  (k-q)t 

i  -  e 


h  '  E  t,  +  ^  -  e 

Fj  =  J  f(z)  dz  +  J  f(z)  dz 

E  tj  +  E 


lim  £  — >  0 


t1+E  h+tj  +  e 

F,=  Jf(z)dz+Jf(z)dz  liimE_i0 

t,  -  E  tj  +  tj  -  E 


f[z)  =  ^le'>,Z 

dz 


2T|/qU(z)  +  ^-| 


t! 


(2.8) 


and  Aac(Nt)  is  the  amplitude  of  the  acoustic  wave  entering  the  substrate,  U(z)  is  the  displacement 
velocity  which  has  the  Bloch  form  U(z)  =  U(z+nt),  E(z)  is  the  electric  field  which  also  is  of  the 
Bloch  form,  e  is  the  piezoelectric  stress  constant,  q  is  the  photon  wavevector,  k  is  the  phonon 
wave  vector,  ti  and  t2  are  the  thickness  of  the  respective  layers,  and  t  =  ti  + 12.  The  M  coefficient 
in  Eq.  (2.7)  is  a  maximum  when  ksq  and  takes  on  the  value  N.  The  Ft  term  represents  the 
contribution  from  the  bulk  and  is  usually  negligible  because  for  Xei  » t,  E(z)  =  const,  and  e(z)  is 
constant  within  each  layer.  The  main  contribution  comes  from  F2  which  represents  the 
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contribution  from  the  discontinuities  at  the  boundaries.  If  the  change  in  permittivity  between  the 
layers  is  small,  E(z)  is  approximately  constant,  and  the  differential  term  in  f(z)  can  be  approximated 
by  Ae.  A  more  complete  expression  for  Eq.  (2.6)  is: 


P  °c 

a c 


(2.9) 


indicating  that  both  the  number  of  layers  and  the  difference  in  the  piezoelectric  constant  between  the 
layers  affect  the  output  power. 

The  reason  for  the  high  conversion  efficiencies  in  a  multilayer  transducer  can  be  seen  by 
looking  at  the  electromagnetic  and  acoustic  dispersion  curves  for  a  multilayered  media.  When  q  = 
k,  the  acoustic  and  electromagnetic  wave  are  phase  matched,  this  phase  matching  occurs  whenever, 
in  the  reduced  zone  scheme,  the  photon  dispersion  curve  intersects  the  phonon  dispersion  curve. 

At  these  points,  the  direct  conversion  of  photons  to  phonons  is  permitted.  The  acoustic  dispersion 
relation  of  a  multilayered  media  in  the  continuum  limit  is  given  by:14 


cos  (qt)  =  cos 


1  +P 


a e. 


2p 


sin 


ac 


(2.10) 


where  p^  =  viP/v2P2>  vi  and  V2  are  the  sound  velocities  and  pi  and  p2  are  the  densities 
of  their  respective  layers.  In  the  continuum  limit,  the  electromagnetic  dispersion  relation  of 
a  multilayered  media  is  given  by:14 


cos  (kt)  =  cos 


1+p; 

2Pel 


sin 


sin 


(2.11) 


where  pel  =  zaiJzal.,  £i  and  £2  are  the  permittivities,  pi  and  \ii  are  the  permeabilities  and  ci  and 
C2  are  the  propagation  velocities  of  electromagnetic  radiation  in  the  respective  medium. 
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Since  the  layer  thicknesses  in  a  millimeter- wave  ZnO  multilayer  transducer  are  on  the  order  of 
a  thousand  Angstroms  or  less,  sin  =  o)  =  0  and  cos  (^=  o)  =  1,  the  electromagnetic  dispersion 

is  a  near  vertical  line  with  k  =  0  for  all  frequencies  in  the  millimeter-wave  region.  If  the  layers  have 
similar  acoustic  properties,  Pac  =»  l,the  acoustic  dispersion  relation  becomes: 


cos  (qt)  =  cos 


(2.12) 


Figure  2.5  is  a  plot  of  the  electromagnetic  and  acoustic  dispersion  curves  for  an  ideal  multilayer 
transducer  with  a  center  frequency  of  32  GHz.  Because  of  the  zone  folding  introduced  by  the 
periodicity  of  the  multilayer  transducer,  it  should  also  be  possible  to  generate  acoustics  at  all  the 
harmonics  of  the  fundamental  frequency  through  umklapp  processes.  Because  the  F2  term  in  Eq. 
(2.7)  is  the  major  source  of  the  acoustic  waves,  gain  occurs  only  at  the  odd  harmonics  of  the 
fundamental  as  can  be  seen  in  Fig.  2.6.  Fig.  2.6(a)  is  the  fundamental  resonance,  the  waves 
generated  at  each  boundary  are  in  phase  and  constructively  interfere  to  produce  the  amplified 
acoustic  wave.  The  second  harmonic,  Fig.  2.6(b),  suffers  destructive  interference  at  each 
boundary  and  the  resulting  wave  is  negligibly  small.  The  third  harmonic,  Fig.  2.6(c),  is  once 
again  in  phase  with  the  wave  generated  in  the  next  interface  and  the  constructive  interference 
produces  an  amplified  acoustic  wave. 

In  a  non- ideal  transducer,  the  values  for  pac  and  pei  will  not  be  exactly  unity,  resulting  in 
imaginary  wave  vectors  in  certain  regions  of  frequency  space,  Fig.  2.7.  These  frequencies  will 
suffer  high  attenuation  rates  and  the  conversion  efficiency  in  these  regions  will  be  quite  small. 

Fig.  2.7  also  indicates  that  even  in  a  non-ideal  multilayer  transducer  there  is  still  direct 
photon-to-phonon  conversion  at  the  intersection  of  the  dispersion  curves  but,  that  in  the  non-ideal 
multilayer  transducer  there  will  also  be  frequencies  where  the  conversion  efficiency  will  be 
negligibly  small.  Non-ideal  multilayer  transducers  should  still  have  the  same  higher  harmonic 
performance  as  the  ideal  transducer. 
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Figure  2.5  Acoustic  and  electromagnetic  dispersion  curves  for  an  ideal  multilayer 
transducer. 
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Figure  2.6  Interference  eff  ects  at  higher  harmonics  in  a  multilayer  transducer. 

(a)  Amplification  at  fundamental  due  to  constructive  interference,  (b)  No  amplification  at 
the  second  harmonic  due  to  destructive  interference,  (c)  Amplification  at  third  harmonic 
due  to  constructive  interference. 
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Figure  2.7  Acoustic  and  electromagnetic  dispersion  curves  for  a  non-ideal  multilayer 
transducer. 
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3.  Results 

The  theoretical  models  of  the  multilayer  transducer  usually  contain  an  identical  number  of 
piezoelectrically  active  and  inactive  layers.  In  practice,  an  odd  number  of  layers  is  used  and  an 
active  layer  is  used  for  both  the  first  and  last  layers.  The  reason  for  using  an  active  layer  as  both 
the  first  and  last  layer  is  that,  as  a  detailed  analysis  given  later  in  the  chapter  will  show,  such  a 
structure  offers  a  slight  improvement  in  the  conversion  efficiency. 

The  optimal  performance  of  the  transducer  is  obtained  when  the  top  electrode  introduces  no 
acoustic  loading  at  the  top  layer  of  the  ZnO.1  The  performance  can  be  optimized  by  making  the  top 
electrode  thickness  much  less  than  an  acoustic  wavelength  or  by  making  the  thickness  a  multiple  of 
one-half  acoustic  wavelength.  The  combination  of  the  high  frequencies  and  slow  sound  speed  of 
the  gold  top  electrode  make  it  practically  impossible  to  make  electrodes  that  are  substantially  thinner 
than  an  acoustic  wavelength.  Using  a  high  multiple  of  the  acoustic  wavelength  is  undesirable 
because  the  top  electrode  resonance  will  narrow  the  bandwidth  of  the  transducer.  A  good 
compromise  for  the  top  electrode  thickness  appears  to  be  between  1000  and  2000A.  Such 
thickness  provided  only  small  top  electrode  loading  without  reducing  the  bandwidth  of  the 
transducer  to  unacceptable  levels.  Gold  is  usually  used  as  the  top  and  counter  electrodes  because 
of  its  good  conductivity  and  ease  of  vacuum  deposition.  In  addition  to  the  gold  electrodes,  a  thin 
layer  of  titanium  is  usually  used  as  an  adhesion  layer  between  the  sapphire  and  the  counter 
electrode,  and  a  thin  layer  of  chromium  is  used  as  an  adhesion  layer  between  the  ZnO  and  the  top 
electrode. 

In  most  transducers,  the  top  electrode  is  a  circular  metallic  region  which  can  be  bonded 
directly  with  any  of  the  standard  thermo-compression  techniques.  Because  of  the  expected 
conversion  efficiencies  in  the  millimeter-wave  frequency  range,  the  transducer  radius  is  made  as 
small  as  possible,  only  20  |im,  too  small  for  direct  bonding  of  the  approximately  50  Jim  gold 
wire.  Instead  of  bonding  directly  to  the  top  electrode  of  the  transducer,  an  "offset  top  dot" 
technique  is  used.15  In  the  offset  top  dot  technique,  a  hole  is  etched  in  the  counter  electrode  before 
the  ZnO  is  deposited.  The  top  electrode  is  deposited  as  usual,  but  then  instead  of  the  circular  top 
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dot  normally  used,  a  specially  shaped  top  dot  is  deposited,  Fig.  2.8.  If  the  large  bond  pad  is 
situated  over  the  region  without  a  counter  electrode,  it  will  not  introduce  significant  capacitive 
loading.  The  transducer,  on  the  other  hand,  is  positioned  over  the  counter  electrode  which 
completes  the  circuit.  With  the  offset  top  dot  technique,  transducers  with  diameters  smaller  than 
the  diameter  of  the  gold  ball  bond  can  be  used. 

The  transducers  are  mounted  in  a  holder  and  the  top  electrode  is  connected  to  the  microstrip 
line  with  a  characteristic  impedance  of  50  Q  by  pressing  the  50-|im-diameter  gold  wire  into  some 
indium  solder.  The  length  of  the  gold  wire  is  kept  very  short  (typically  less  than  0.5  mm)  to 
minimize  the  series  reactance.  In  the  Ka-Band,  a  coaxial  connector  at  the  other  end  of  the  holder  is 
connected  to  the  microstrip  line  and  is  used  to  bring  the  signals  into  and  out  of  the  transducer.  Fig. 
2.9(a).  Because  of  the  high  loss  encountered  using  coaxial  connectors  in  the  W-band,  a  transition 
between  the  WR-10  wave  guide  (characteristic  impedance  400  Q )  and  the  40  Q  microstrip  line 
was  constructed  using  a  tapered  reduced-height  ridge-waveguide  section16,  Fig.  2.9(b).  The 
conversion  efficiency  of  the  transducer  is  measured  using  a  pulse  echo  technique  in  which  a  short 
rf-tone  burst  is  applied  to  the  transducer  and  the  resulting  echo  off  the  opposite  end  of  the  sapphire 
flat  are  measured.  In  order  to  improve  the  dynamic  range  of  the  measurements,  a  superheterodyne 
detection  scheme  in  which  the  return  signal  is  mixed  down  and  detected  at  a  lower  frequency  is 
employed. 

(1)  Verification  of  Theoretical  Predictions 

In  the  beginning  section,  it  was  stated  that  for  most  transducers  it  is  the  electrical  response 
that  dominates  the  conversion  loss  of  the  transducer.  The  validity  of  this  assumption  in  this 
frequency  range  can  be  tested  by  changing  the  transducer  configuration  and  measuring  the 
conversion  efficiency.  If  the  conversion  loss  varies  according  to  Eq.  (2.2),  the  electrical  losses  are 
the  dominant  factor,  if,  on  the  other  hand,  the  response  deviates  from  Eq.  (2.2)  the  acoustic  losses 
must  be  considered  in  evaluating  the  multilayer  transducer.  The  easiest  way  to  modify  the  electrical 
conversion  loss  is  to  alter  the  real  part  of  the  transducer  impedance,  Ra.  Because  the  capacitance  of 
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bond  pad 


Figure  2.8  Schematic  of  bonding  to  the  transducer  with  the  offset  top  dot  technique, 
(a)  top  view,  (b)  side  view. 
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Tuning  stub 


Figure  2 .9  Electrical  connections  to  the  transducers, 
(a)  Ka-band,  (b)  W-band. 
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the  transducer  is  proportional  to  the  area  of  the  transducer,  Ra  is  most  dramatically  affected  by 
changing  the  radius  of  the  transducer. 

In  order  to  determine  the  effect  of  transducer  radius  on  conversion  efficiency,  two 
transducers,  both  with  offset  top  dots,  were  deposited  on  the  same  sapphire  flat  using  the 
parameters  described  in  Table  2.1. 17  The  first  transducer  had  a  radius  of  20  fim  while  the  second 


Total  number  of  layers 

15 

ZnO  layer  thickness 

< 

o 

in 

On 

Top  electrode  thickness 

1500  A  Ti/Au 

Counter  electrode  thickness 

1000  A  Cr/Au 

Transducer  radius 

20  pm 

Sapphire  substrate  thickness 

2.6  mm 

Table  2.1.  Transducer  parameters  for  Ka-Band  transducers. 


transducer  had  a  radius  of  60  pm.  The  highest  conversion  efficiency  for  both  transducers  occurred 
near  28  GHz  with  the  20  pm  radius  transducer  having  a  "raw"  conversion  efficiency  of  18  dB 
higher.  The  diffraction  losses  can  be  calculated  to  be  3.1  dB  for  the  20  pm  radius  transducer  and 
1. 1  dB  for  the  60  pm  radius  transducer  at  28  GHz.  When  the  diffraction  losses  are  subtracted 
from  the  data,  the  difference  between  the  two  transducers  is  20  dB.  Other  minima  can  also  be 
correlated  between  the  two  transducers.  The  minima  near  30  and  35  GHz  show  a  diffraction 
corrected  conversion  efficiency  difference  of  18  dB  as  shown  in  Fig.  2.10.  If  the  conversion  loss 
is  calculated  using  Eqs.  (2.2)  and  (2.4),  the  expected  difference  between  the  two-way  conversion 
loss  is  19  dB.  This  close  agreement  between  theoretical  and  experimental  values  indicates  that  it 
is  the  losses  in  the  electrical  system,  rather  than  the  losses  in  the  acoustic  system,  that  dominate  the 
untuned  conversion  efficiency. 

One  of  the  conclusions  that  can  be  drawn  from  Eqs.  (2.12)  and  (2.13)  is  that  the  transducers 
should  operate  at  all  the  odd  harmonics  of  the  fundamental  frequency.  To  test  this  prediction,  one 
of  the  Ka-band  transducers  was  operated  in  the  the  W-band  to  see  if  the  transducer  had  measurable 
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Figure  2 .10  Comparison  of  transducer  radius  on  conversion  efficiency, 
(a)  20  Jim  radius  transducer,  (b)  60  Jim  radius  transducer. 
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conversion  efficiency  at  the  third  harmonic.  Fig.  2.1 1  is  a  plot  of  the  conversion  efficiency  in  the 
W-band  and,  for  comparison,  a  plot  of  the  Ka-band  response  multiplied  by  a  factor  of  three  in 
frequency.  The  conversion  efficiency  at  three  times  the  fundamental  Ka-band  resonance  is  useful 
only  in  comparing  relative  responses  as  a  function  of  frequency.  As  can  be  seen  in  the  figure,  the 
response  at  28.6  GHz  in  the  Ka-band  is  reproduced  nicely  by  its  third  harmonic  counterpart  near 
87  GHz. 

Earlier  work  on  single  layer  tilted  ZnO  films  has  shown  that  tilted  ZnO  transducers  still  had 
longitudinal  wave  coupling  factors  that  were  about  half  the  magnitude  of  the  normal  ZnO  value.11 
If  each  of  the  layers  in  the  multilayer  transducer  has  some  longitudinal  coupling,  there  will  be  some 
"bulk"  response  from  the  transducer  near  the  half-wavelength  resonance  frequency  for  the  whole 
transducer.  Measurements  were  carried  out  on  two  of  the  Ka-band  transducers  in  the  5-12  GHz 
frequency  range  to  determine  the  bulk  response.  The  conversion  loss  for  the  transducer  mounted 
in  hole  number  4  is  plotted  in  Fig.  2.12.  The  highest  conversion  efficiency  is  achieved  around 
9  GHz  and  is  probably  the  second  harmonic  of  the  fundamental  near  4  GHz.  The  bulk  response 
of  the  transducer  indicates  that  the  value  of  kbulk  =  0.21. 

Because  the  offset  top  dot  technique  had  never  been  used  at  millimeter-wave  frequencies 
before,  its  effectiveness  with  such  short  wavelength  radiation  needed  to  be  evaluated.  Two 
transducers  were  fabricated  on  the  same  substrate  using  the  parameters  in  Table  2.1  except  that  the 
transducer  radius  was  increased  to  60  |im.  One  transducer  had  a  circular  top  dot  with  the  gold  ball 
bonded  directly  to  the  top  dot,  the  other  had  an  offset  top  dot.17  The  two  transducers  had  similar 
conversion  efficiencies  and  frequency  responses  as  can  be  seen  in  Fig.  2. 13.  This  result  suggests 
that  the  loading  of  the  transducer  by  the  gold  ball  and  wire  is  probably  a  secondary  effect  and  does 
not  change  the  response  of  the  transducer.  It  also  appears  that  the  effects  of  the  unusual  electrical 
connection  used  in  the  offset  top  dot  technique  are  negligible  at  these  frequencies  even  when  the 
electrical  wavelengths  are  only  a  few  millimeters. 
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Fig  2 .1 1  Frequency  response  of  a  multilayer  transducer  operated  at  the  third  harmonic. 
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Figure  2 .12  Bulk  response  of  multilayer  transducer. 
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Figure  2 .13  Comparison  of  conversion  efficiencies  of  the  two  top  electrode  bonding 
techniques. 

(a)  offset  top  electrode,  (b)  direcdy  bonded  top  electrode. 
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4.  Piezoelectric  Modulation 

The  theoretical  foundation  for  the  multilayer  transducer  requires  that  there  exist  some  periodic 
modulation  of  the  piezoelectric  constant  in  order  to  achieve  gain  over  the  single  layer  transducer. 
The  foundation  for  the  multilayer  ZnO  transducer  postulated  that  the  normal  and  tilted  layers 
would  have  a  different  orientation  of  the  c-axis  and  thus  a  different  piezoelectric  coupling  constant 
for  longitudinal  waves.  Early  work  on  films  that  were  several  microns  thick  had  shown  that  the 
c-axis  of  the  bulk  crystal  could  be  rotated  by  depositing  the  ZnO  at  a  given  angle.13  It  has  been 
shown  that  the  change  in  c-axis  orientation  does  not  occur  abruptly  but  gradually  changes  over 
some  finite  transition  distance.18  If  the  distance  required  for  the  c-axis  to  change  orientation  is  on 
the  order  of  an  acoustic  wavelength,  the  theoretical  approach  used  in  Eq.  (2.12)  may  not  accurately 
model  the  physical  situation.  Additionally,  an  interface  layer  might  necessitate  changes  in  the 
thickness  of  either  the  normal  or  tilted  layers  in  order  to  get  the  correct  periodicity  for  maximum 
acoustic  power  generation.  Experimental  work  on  thin  film  tilted  ZnO  found  that  the  films  still  had 
substantial  longitudinal  wave  coupling  even  though  the  entire  film  was  deposited  at  40°. 11  The 
authors  suggested  that  it  was  possible  that,  rather  than  having  a  "tilted"  c-axis  orientation  as 
predicted,  the  film  had  more  of  an  amorphous  character.  In  order  to  more  closely  examine  the 
grain  structure  of  the  layers  and  transition  distance  between  crystal  structure,  the  transducers  were 
examined  using  electron  microscopy.19 

The  transducers  were  first  examined  with  the  scanning  electron  microscope  (SEM).  The 
sapphire  substrates  were  first  cut  from  the  back  side  to  within  0.005  inches  of  the  top  surface.  The 
remaining  sapphire  and  the  transducer  were  then  cleaved  by  breaking  the  substrate  along  the  cuts. 
This  technique  has  proven  effective  in  examining  single  layer  transducers  and  seems  to  leave  the 
exposed  cross  section  free  of  debris  and  damage.  The  samples  are  then  mounted  and  sputter 
coated  with  100-200  A  of  gold.  The  SEM  micrographs  of  the  multilayer  structures  reveal  that 
every  other  layer  of  the  transducer  appears  bright  Fig.  2.14  (micrograph  courtesy  of  Rosemarie 
Koch,  Center  for  Materials  Research,  Stanford  University).  One  possible  explanation  of  the  SEM 
micrographs  is  that  the  "normal"  and  "tilted"  layers  cleave  at  different  angles  yielding  a  "stair  step" 
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Figure  2. 14  SEM  micrograph  of  a  multilayer  transducer. 
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topography  that  is  imaged  as  bands  by  the  SEM.  In  some  high  magnification  images,  the  transition 
between  the  normal  and  tilted  layer  appears  rather  gradual  and  the  transition  between  the  tilted  and 
normal  layer  seems  sharper.  Such  transitions  suggest  that  the  change  in  the  crystal  axis  orientation 
from  a  tilted  film  to  a  normal  film  may  occur  over  a  very  short  distance,  but  that  the  change  in  the 
c-axis  orientation  from  normal  to  tilted  layers  occurs  over  a  significantly  longer  distance.  The  SEM 
is  unable  to  precisely  determine  the  behavior  of  the  interface  between  the  layers  because  the  sample 
must  be  coated  with  at  least  100A  of  Au  to  provide  enough  contrast 

In  order  to  observe  the  interface  between  the  layers  without  the  gold  coating,  samples  were 
prepared  for  examination  using  the  transmission  electron  microscope  (TEM).  For  the  SEM 
samples,  actual  transducers  whose  conversion  loss  had  been  measured  were  imaged. 
Unfortunately,  TEM  thin  film  cross  sections  could  not  be  made  from  films  deposited  on  sapphire 
substrates  because  of  the  extreme  hardness  of  sapphire.  Silicon  was  chosen  as  the  substrate  for  the 
TEM  samples  because  other  work  had  shown  that  high  quality  ZnO  could  be  grown  on  silicon20 
and  silicon  was  soft  enough  for  use  in  TEM  thin  sections.  Unfortunately,  commerically  available 
silicon  wafers  are  so  thin  that  acoustic  measurements  could  not  be  performed  on  the  transducers, 
so  the  conversion  efficiencies  of  transducers  examined  with  the  TEM  are  unknown.  The 
transducer  sample  was  deposited  on  a  Si  wafer  with  a  1000A  thermal  oxide.  Because  the  sample 
did  not  need  to  resonate  in  a  particular  frequency  band,  convenient  parameters  were  chosen  for  the 
deposition  thicknesses  as  shown  in  Table  2.2. 


Total  number  of  layers 

10 

ZnO  layer  thickness 

iooo  A 

Top  electrode 

none 

Counter  electrode  thickness 

1000  A  Cr/Au 

Substrate 

Si  with  1000  A 

thermal  Si02 

Table  2.2  Parameters  for  the  transducer  examined  by  the  TEM. 

The  TEM  micrographs  indicate  the  same  banding  structure  as  the  SEM  images  but  the 
individual  layers  can  clearly  be  distinguished  as  seen  in  Fig.  2.15  (micrograph  courtesy  of  Ann 
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Figure  2.15  TEM  micrograph  of  a  multilayer  transducer. 
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Marshall,  Center  for  Materials  Research,  Stanford  University).  The  banding  is  indicative  of 
crystal  changes  between  the  normal  and  tilted  layers.  The  TEM  micrographs  indicate  that  the 
"normal"  layer  is  1050A  thick,  near  the  expected  value  of  lOOOA.  The  "tilted"  layers  are  750A 
thick,  significantly  less  than  the  expected  value  of  1000A.  The  error  in  the  TEM  measurement  is 
less  than  ±5%.  The  discrepancy  in  the  thickness  of  the  tilted  layer  could  come  from  several 
factors.  The  tilted  layer  may  initially  grow  more  slowly  than  the  rate  calculated  using  thicker  film 
would  indicate.  It  is  also  possible  that  some  of  the  thickness  attributed  to  the  normal  layer  is 
deposited  during  the  tilted  deposition  phase  but  its  crystal  structure  more  closely  resembles  the 
normal  layer  and  is  imaged  as  such  by  the  TEM. 

The  TEM  micrographs  do  not  show  the  large  scale  changes  in  crystal  grain  structure  between 
the  normal  and  tilted  layers  that  has  been  seen  with  thicker  films.9  Electron  diffraction  analysis  of 
the  TEM  samples  indicates  that  both  the  normal  and  tilted  layers  show  strong  c-axis  orientation 
parallel  to  the  surface  normal  indicating  that  the  distance  required  for  the  crystal  axes  to  change 
orientation  is  at  least  750A.  The  diffraction  analysis  also  shows  no  significant  differences  between 
layers  that  were  deposited  at  ±  40°  to  the  surface  normal.  The  small  change  in  the  c-axis  orientation 
through  the  entire  multilayer  transducer  indicates  the  c-axis  orientation  is  not  the  major  cause  of  the 
piezoelectric  modulation.  Even  though  the  change  in  crystal  orientation  between  the  two  layers 
appears  small,  the  sharpness  of  the  transition  is  still  important  because  that  will  set  an  upper  limit 
on  the  frequencies  that  can  be  generated  using  this  method.  Detailed  examination  of  the  interface 
between  the  normal  and  tilted  layers  indicates  that  the  change  in  the  TEM  image  contrast  occurs 
over  a  distance  of  approximately  50  A. 

Another  mechanism  that  could  cause  changes  in  the  piezoelectric  constant  is  differences  in  the 
chemical  composition  between  the  layers.  A  transducer  with  the  same  parameters  as  those  in  Table 
2.2,  with  the  substrate  changed  to  sapphire,  was  examined  using  Auger  electron  spectroscopy.  A 
depth  profile  through  transducer  is  constructed  by  using  argon  ions  to  etch  a  hole  in  the  sample 
while  the  Auger  electrons  from  the  region  of  interest  are  monitored.  Using  this  method,  a  profile  of 
elemental  composition  versus  depth  is  created  as  shown  in  Fig.  2.16.  The  electron  count  rate  can 
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Figure  2 .16  Depth  profile  of  a  multilayer  transducer  using  Auger  electron  analysis. 
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then  be  calibrated  against  a  sample  of  single  crystal  ZnO  to  determine  if  the  layers  vary  in  chemical 
composition.  The  spectroscopy  reveals  that  there  are  no  compositional  changes  that  occur  at  the 
layer  boundaries.  The  Auger  analysis  did  yield  one  anomalous  result;  the  first  2000A  of  ZnO 
appears  to  be  deposited  with  a  higher  oxygen  concentration  than  the  layers  that  are  deposited 
afterward.  The  later  layers  have  the  same  stoichiometry  as  single  crystal  ZnO.  The  abundance  of 
oxygen  present  in  the  beginning  of  the  deposition  probably  results  from  oxygen  adhering  to  the 
surface  of  the  Zinc  target  or  excess  oxygen  adhering  to  the  surface  of  the  counter  electrode.  The 
higher  oxygen  concentration  may  degrade  the  performance  of  the  transducer  but  should  not 
interfere  with  the  piezoelectric  modulation. 

5.  Computer  Simulations 

The  effects  of  various  physical  parameters  on  transducer  performance  have  been  investigated 
theoretically  using  a  computer  simulation  based  on  the  earlier  work  of  Hadimioglu3  and  Sittig.20 
The  multilayer  transducer  can  be  modeled  by  a  Mason's  equivalent  circuit  by  making  the  following 
assumptions: 

1)  The  transducer  is  considered  to  be  one  dimensional  in  the  direction  of  sound  propagation. 

The  transducer  radius  is  much  larger  than  an  acoustic  wavelength  and  the  electric  and  acoustic 
fields  vary  only  in  the  thickness  direction. 

2)  Each  of  the  layers  may  have  acoustic  losses  that  vary  as  a  function  of  frequency,  but  the 
losses  in  each  layer  must  vary  in  an  identical  manner.  The  acoustic  losses,  a,  must  be  small 
such  that  a  «  w/c 

3)  Each  layer  is  considered  electrically  lossless. 

4)  The  top  surface  is  approximated  as  a  free  surface  and  the  substrate  surface  is  approximated  as 
a  real  impedance  load. 

5)  Only  the  longitudinal  acoustic  modes  generated  in  each  layer  are  considered. 

With  the  above  assumptions,  a  series-connected  multilayer  transducer  can  be  modeled 
by  the  Mason  circuit  in  Fig.  2.17.20  The  following  definitions  are  used  in  the  Mason's 
model: 


34 


CELLO  CELLO  ....  CELL  N 

LAYER  0  LAYER  1  LAYER  M 

TOP  w  r  »  COUNTER 

ELECTRODE  v00-  r00-  *00  v01- r01-  *01  VNM- rNM- «NM  ELECTRODE 
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Fig  2 . 17 Mason's  equivalent  circuit  for  a  multilayer  transducer. 


v  is  the  longitudinal  sound  velocity, 

r  is  the  mechanical  impedance, 

t  is  the  layer  thickness. 


0  =fi-+ya(f)t 

<t>  =  0)^33 


is  the  clamped  static  capacitance  of  the  transducer,  where  A  is  the 
area  of  the  transducer  and  e  is  the  dielectric  constant  of  the  medium, 
is  the  propagation  length  in  the  medium,  where  a(f)  is  the  loss  in  the 
medium  which  can  be  a  function  of  frequency, 
is  the  coupling  constant  of  the  medium,  where 


hoo  =  ^21  =  k  /£2l]  '2  is  the  piezoelectric  coupling  constant,  633  is  the 

e§3  \ey 

piezoelectric  stress  constant,  is  the  permittivity  at  constant  strain,  C33 

is  the  stiffness,  and  kt  is  the  electro-mechanical  coupling  constant, 
the  subscripts  used  in  the  Mason's  model  have  the  following  meaning: 
te  top  electrode, 
ce  counter  electrode, 
sub  substrate, 
ij  layer  j  within  cell  i. 

In  the  Mason's  model,  the  current  I,  the  voltage  V,  the  force  F  and  the  displacement  velocity 
U  of  one  layer  are  related  to  the  same  quantities  of  the  next  layer  by  the  following  relation: 


rF  ,  1 

F  ] 

n+l 

n 

U„+l 

=(\] 

u 

n 

.V„+l. 

V 

L  n  J 

(2.13) 
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Because  the  layers  are  electrically  excited  in  series  configuration,  the  current  I  is  the  same  for  all 
layers.  The  components  of  the  A  matrix  can  be  derived  from  elementary  circuit  analysis  of  the 
Mason's  model  in  Fig.  2.17. 
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j  sin  0 


hsin  0 
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G)Cn 
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(2.14) 


where  co  is  the  angular  frequency,  Z  is  the  impedance,  Co  is  the  clamped  capacitance,  and  h  the 
piezoelectric  coupling  constant.  Each  layer  will  have  a  different  A  matrix  that  depends  on  the 
material  parameters  for  that  layer.  The  A  matrices  for  the  top  and  counter  electrodes  can  be  created 
by  setting  their  piezoelectric  coupling  constant,  h,  to  zero. 

Eq.  (2.13)  can  be  broken  into  parts  so  that  the  force  F  and  the  displacement  velocity  U  can  be 
solved  iteratively. 
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and  thus: 
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where  Ay  refer  to  the  elements  of  the  A  matrix  in  Eq.  (2.14).  In  this  way,  the  forces  and  velocities 
at  the  top  surface  can  be  related  to  the  forces  and  velocities  at  the  counter  electrode-substrate 
interface.  The  boundary  conditions  at  these  interfaces  yield  a  set  of  simultaneous  equations  that  are 
solved  for  the  acoustic  power  radiated  into  the  substrate. 
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Typical  physical  parameters  used  as  input  for  the  computer  simulation  are  shown  in  Table 


2.3. 


Normal  layer  velocity 

5800  m/s 

Tilted  layer  velocity 

5500  m/s 

Normal  layer  coupling  constant 

.2$ 

Tilted  layer  coupling  constant 

.21 

Normal  layer  impedance 

3.6  x  lO'7 

Tilted  layer  impedance 

3.2  x  10-7 

Top  electrode  velocity 

2800  m/s 

Counter  electrode  velocity 

3100  m/s 

Table  2.3.  Transducers  parameters  used  in  computer  simulation. 


The  values  in  Table  2.3  are  derived  from  experimental  and  theoretical  work  done  on  single  layer 
transducers  at  lower  frequencies  and  provide  a  starting  point  for  estimating  the  physical  parameters 
at  higher  frequencies.  In  the  simulation,  the  values  of  the  various  physical  parameters  and  layer 
thicknesses  are  changed  until  a  reasonable  fit  to  the  experimental  data  is  obtained.  Unfortunately, 
computer  simulation  using  this  model  has  yielded  only  satisfactory  agreement  with  the 
experimental  results,  Fig.  2. 18.  The  simulation  does  not  consider  several  points  that  may  have  a 
bearing  on  the  measured  value.  Because  the  program  considers  only  the  longitudinal  wave 
propagation,  the  effects  of  any  shear  waves  that  may  be  generated  in  the  "tilted"  layers  and  the 
effects  they  have  on  the  longitudinal  waves  generated  in  the  "normal"  layers  is  not  addressed  by  the 
simulation. 

Experimentally,  no  evidence  of  shear  wave  generation  was  found  at  powers  20  dB  lower 
than  the  two-way  longitudinal  wave  power  in  either  the  Ka  or  W-Band.  If  the  effects  of  the  shear 
waves  are  small,  it  should  be  possible  to  model  the  shear  wave  effects  by  modifying  the  velocity, 
impedance,  and/or  coupling  constant  of  the  tilted  layer,  but  reasonable  fits  were  obtained  only 
when  the  estimated  physical  parameters  differed  greatly  from  their  bulk  values.  It  may  be  possible 
that  while  the  shear  wave  power  entering  the  substrate  is  small,  shear  waves  play  an  important  role 
in  the  longitudinal  wave  generation. 
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The  simulation  also  assumes  that  each  layer  is  uniform  in  crystal  orientation  and  is  identical  to 
all  the  other  layers  of  that  orientation.  While  all  the  layers  are  deposited  under  nearly  identical 
conditions,  it  is  possible  that  the  first  normal  layer  which  is  deposited  on  gold  could  have  different 
physical  properties  than  the  other  normal  layers  which  are  deposited  on  tilted  ZnO.  The  use  of 
both  positions  II  and  in  for  depositing  the  tilted  ZnO  may  result  in  some  differences  between  the 
tilted  layers.  Several  simulations  were  run  with  different  physical  parameters  for  the  ZnO  layers 
deposited  in  position  II  and  in,  but  good  agreement  with  the  experimental  data  required  highly 
atypical  values. 

Another  aspect  of  the  multilayer  transducers  that  has  been  difficult  to  model  is  the  transition 
region  between  the  normal  and  tilted  layers.  The  TEM  analysis  indicated  that  the  distance  between 
the  ZnO  crystal  structures  occurred  over  a  short  distance,  but  the  interface  layer  could  still  play  an 
important  role  in  the  longitudinal  wave  generation.  Two  methods  have  been  used  to  simulate  a 
small  transition  region  between  the  layers  on  the  computer.  The  first  method  assigned  the 
transition  region  to  one  or  the  other  of  the  layers  and  used  the  physical  properties  of  the  assigned 
layer  for  the  transition  layer.  The  second  method  involved  modeling  the  reorientation  region  as  a 
thin  layer  between  each  of  the  two  larger  layers.  The  physical  parameters  of  the  transition  region 
were  computed  to  be  the  arithmetic  mean  of  the  layers  on  either  side.  Neither  approximation 
yielded  significantly  better  agreement  with  experimental  data  than  simply  omitting  the  transition 
region  altogether. 

6.  Extension  to  Higher  Frequencies 

In  the  past,  coherent  millimeter-wave  frequency  acoustic  wave  generation  has  been 
accomplished  at  frequencies  as  high  as  1 14  GHz  by  placing  a  quartz  crystal  in  the  high  electric  field 
region  of  a  cavity.21  Because  the  W-band  multilayer  ZnO  transducers  have  a  conversion  loss  that 
is  typically  50  dB  less  than  the  quartz  transducers,  it  is  useful  to  examine  the  possibilities  of 
operating  the  ZnO  multilayer  transducers  at  even  higher  frequencies.  Several  practical  problems 
must  be  overcome  if  the  higher  frequency  transducers  are  to  be  constructed  and  evaluated.  The 
acoustic  attenuation  in  sapphire  at  higher  frequencies  must  be  considered,  and  while  no  one  has 
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measured  the  acoustic  attenuation  of  sapphire  millimeter- wave  frequencies  higher  than  100  GHz, 
results  on  quartz  indicate  that  unless  the  sample  is  cooled  to  below  20  K  the  acoustic  attenuation 
may  pose  a  serious  problem.21  The  W-band  multilayer  transducers  have  been  operated  at 
temperatures  as  low  as  2  K  without  any  damage  to  the  transducer  indicating  that  the  transducers 
can  be  cooled  to  very  low  temperatures  to  reduce  the  attenuation  if  necessary.  Cooling  the  sample 
should  also  reduce  the  acoustic  attenuation  in  the  ZnO,  leading  to  higher  conversion  efficiencies. 

Another  major  problem  is  that  the  current  pulse  echo  detection  technique  places  extremely 
tight  specifications  on  the  parallelism  and  flatness  of  the  sapphire  substrates.  Any  surface 
irregularities  or  imperfections  on  the  substrates  will  also  place  an  upper  limit  on  the  operating 
frequency  that  can  be  detected  using  this  technique.  In  the  W-band,  the  non-flatness  and 
non-parallelism  of  the  sapphire  substrate  were  already  affecting  the  measurements.  At  300  GHz, 
the  wavelength  of  sound  in  sapphire  will  be  about  350  A  and  would  require  extremely  high 
tolerances  on  the  sapphire.  Even  if  the  transducer  radius  is  as  small  as  20  pm,  the  sapphire 
substrate  would  still  have  to  be  flat  to  within  3  pm  over  the  1/4"  diameter  of  the  sapphire  rod  to 
ensure  destructive  interference  did  not  completely  destroy  the  incoming  signal. 

One  way  to  overcome  the  substrate  problems  is  to  detect  the  signal  with  some  type  of 
incoherent  detector  which  would  not  be  as  sensitive  to  the  substrate  condition.  Other  researchers 
have  shown  that  if  the  sound  can  be  coupled  from  the  substrate  into  superfluid  helium  under  the 
right  conditions,  the  helium  is  capable  of  propagating  the  200  GHz  sound  for  distance  of 
millimeters.22  Using  this  approach,  the  detector  could  be  located  at  some  distance  to  decrease  any 
rf  coupling  between  the  transducer  and  detector. 

Several  problems  also  confront  the  transducer  design  if  the  frequency  range  is  to  be 
increased.  One  of  the  problems  that  may  be  encountered  in  extending  the  operating  frequencies 
above  100  GHz  is  maintaining  the  periodicity  in  the  piezoelectric  modulation.  In  a  non-ideal 
transducer,  the  layer  thicknesses  are  going  to  be  distributed  around  one-half  wavelength  because  of 
variations  in  the  deposition  parameters.  One  can  estimate  the  effect  of  thickness  variations  on  the 
output  power.  If  the  separation  between  active  layers  is  not  exactly  one  acoustic  wavelength,  the 
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wave  produced  in  that  layer  will  have  some  phase  offset  relative  to  the  entering  wave.  The  acoustic 
amplitude  at  the  substrate  would  then  be: 
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(2.17) 


where  Ao  is  the  acoustic  amplitude  generated  at  each  interface,  and  <J)(n)  is  the  phase  difference 
between  the  "ideal"  wave  and  the  acoustic  wave  generated  at  the  boundary  between  the  layers.  If 
N  is  large,  Eq.  (2.17)  can  be  approximated  by: 
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\c~wjmcos(m  (2.18) 
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Where  P(<J>)  is  a  weighting  function  describing  the  distribution  of  (J).  If  <J>(n)  is  randomly  distributed 
about  0  such  that  -0  ^  <}>  ^  0, 


0,  <J>  <  — 0 


P(<t>)  =  \  — ,  -0  £  <{>  £  +0 
20 

0,  <}>  >  +0 


Substituting  Eq.  (2.19)  into  Eq.  (2.18)  and  integrating  yields: 


(2.19) 
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sin(0) 
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=  AqN  sinc(0) 


(2.20) 


The  effect  on  the  power  radiated  into  the  substrate  can  be  estimated  by  expanding  sinc2(0)  in 
a  Taylor  series  about  0  =  0. 


42 


(2.21) 


sinc2(0)  =  1  -  —  02  -t- 

If  0  is  small,  then  to  first  order,  the  power  entering  the  substrate  will  be  unaffected  by  the 
thickness  variations.  A  good  estimate  for  0  can  be  found  by  estimating  the  error  in  the  thickness 
caused  by  variations  in  the  deposition  rate  and  time  for  each  layer.  The  error  in  the  deposition  is 
probably  limited  to  a  few  monolayers  of  ZnO  or  about  20A.  For  transducers  operating  in  the  Ka- 
Band,  the  acoustic  wavelength  in  ZnO  is  approximately  2000A  and  0  =  7.2°,  yielding  a  two-way 
decrease  in  power  of  only  0.05  dB.  For  transducers  operating  in  the  W-Band,  the  effect  is  more 
pronounced,  but  still  small.  The  typical  wavelength  in  the  W-band  is  approximately  600  A  and  0  = 
24.0°,  resulting  in  a  two-way  power  loss  of  0.5  dB.  The  two-way  conversion  efficiency  would 
not  be  reduced  by  10  dB  until  the  acoustic  wavelength  reaches  150  A  or  about  350  GHz. 

The  minimum  distance  required  to  modulate  the  piezoelectric  constant  between  the  layers  will 
also  determine  the  maximum  operating  frequency  of  multilayer  ZnO  transducers.  If  the 
piezoelectric  constant  changes  over  the  same  region  as  the  TEM  image  contrast,  about  50  A,  it 
should  be  possible  to  deposit  transducers  with  layers  as  thin  as  100  A,  corresponding  to  operating 
frequencies  of  about  300  GHz,  and  still  achieve  reasonable  conversion  efficiencies.  Because  the 
output  power  increases  as  the  number  of  layers  squared,  it  is  also  possible  that  the  number  of 
layers  could  be  increased  by  an  order  of  magnitude  or  more  to  generate  sufficient  power  at  higher 
frequencies. 

In  conclusion,  multilayer  ZnO  acoustic  transducers  have  been  designed  and  tested  at 
operating  frequencies  from  24  GHz  up  to  100  GHz.  These  transducers  are  by  far  the  most 
efficient  transducers  in  this  frequency  range  ever  constructed.  These  transducers  offer  the 
possibility  of  designing  a  superfluid  helium  acoustic  microscope  with  a  resolution  as  high  as  30A. 
It  also  appears  that  multilayer  ZnO  transducers  can  generate  coherent  phonons  at  frequencies  as 
high  as  300  GHz  with  reasonable  conversion  efficiencies.  Although  some  features  of  the 
transducer's  performance  cannot  be  accounted  for  by  computer  simulations,  the  transducer's 
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performance  closely  matches  the  theoretical  predictions  regarding  harmonic  generation  and 
dependence  on  transducer  radius. 

III.  Propagation  Studies  in  Liquid  Helium 

Liquid  helium  cooled  below  the  Lambda  point  is  a  suitable  media  for  propagating  sound  at 
the  frequencies  of  interest  in  this  report  and  it  is  important  to  increase  our  understanding  of  the 
linear  and  nonlinear  properties  of  this  liquid.  The  non-linear  properties  were  investigated  with  a 
spherically  converging  acoustic  beam.  Since  the  power  densities  are  very  high  at  the  focus,  the 
nonlinear  effects  are  enhanced,  particularly  as  the  frequency  is  increased.  The  use  of  very  low 
temperatures  prevents  scattering  by  thermal  phonons,  but  the  signal-to-noise  ratio  is  limited  by  the 
onset  of  noise-stimulated  nonlinear  scattering.  This  three-phonon  process,  which  becomes 
stronger  as  the  frequency  increases,  produces  an  effective  attenuation  as  it  absorbs  power  from  the 
fundamental  frequency  of  the  signal  beam.  This  "attenuation"  can  be  greatly  decreased,  and  thus 
the  signal-to-noise  ratio  increased,  simply  by  applying  pressure  to  the  helium  bath. 

1.  The  Three-Phonon  Process 

Figure  3.1  depicts  a  three-phonon  scattering  event  in  which  two  phonons  collide  and  form  a 
third  phonon  (the  reverse  process  also  occurs).  The  collision  conserves  energy  and  momentum. 
The  phonon  energy  is  hf,  where  f  is  the  frequency  and  h  is  Planck's  constant;  thus  conservation  of 
energy  requires  fi  +  f2  =  fs-  The  phonon  momentum  is  hk/2jt;  thus  conservation  of  momentum 
requires  kj  +  k2  =  k3.  If  the  dispersion  were  linear,  i.e.,  f=ck/2tt  for  all  frequencies  (where  c  is  a 
constant  speed  of  sound),  then  the  two  conservation  equations  require  that  the  phonons  be 
collinear— all  momenta  in  the  same  direction.  If  the  dispersion  is  upward,  as  shown  in  the  Figure, 
then  the  sound  speed  c  is  increasing  with  frequency.  Thus  the  higher  energy  phonon  fs  will  have 
less  momentum  than  the  dispersionless  case.  The  two  input  phonons  must  then  intersect  at  a  finite 
angle  so  that  some  of  their  momentum  will  be  cancelled  and  the  momentum  can  be  conserved  (the 
drawing  in  the  figure  depicts  this  upward  dispersion  case).  If  the  dispersion  is  downward, 
however,  the  higher  energy  phonon  has  too  much  momentum  (more  than  can  be  obtained  from  the 
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FIGURE  3 . 1  Diagram  showing  the  scattering  of  two  phonons  to  form 
a  third  for  various  dispersion  cases. 
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two  smaller  phonons)  so  that  momentum  cannot  be  conserved.  Downward  dispersion  thus  makes 
the  three-phonon  process  much  less  likely. 

It  must  also  be  pointed  out  that  when  the  dispersion  is  linear,  the  three-phonon  process  that 
dominates  is  not  the  difference  frequency  process,  associated  with  gain,  but  the  sum  process, 
associated  with  the  formation  of  higher  harmonics,  and,  in  normal  liquids,  a  shock  wave.  Studies 
of  harmonic  formation  have  been  made  for  focused  acoustic  beams  in  cryogenic  liquids.23 

The  three-phonon  process  is  frequency  dependent,  with  an  effective  attenuation  that  goes  as 
the  fifth  power  of  frequency.24’25’26  Even  worse,  once  a  scattering  event  occurs  in  which  the 
beam  phonon  is  scattered  into  two  smaller  phonons,  these  products  are  exactly  phase-matched  with 
all  the  other  signal  phonons  to  scatter  them,  in  an  avalanche-type  process,  as  demonstrated  through 
experiments  using  plane  waves.27  In  optics,  the  stimulated  amplification  of  spontaneous  decay 
products  has  long  been  understood.  It  was  first  demonstrated  in  visible  light  by  Akhmanov 
et  a/.,28  Magde  and  Mahr,29  and  Harris  et  a/. 30.3 1,3 2  jn  the  last  experiment,  for  example,  just  by 
shining  a  high  frequency  beam  onto  a  lithium  niobate  crystal  the  spontaneous  decay  products  could 
actually  be  amplified  up  to  where  they  were  visible  with  the  unaided  eye.  In  the  acoustic 
experiments,  Foster  and  Putterman27  showed  that  the  process  could  go  on  in  the  higher 
nonlinearity  of  the  helium  until  the  input  beam  ("pump")  is  actually  depleted  by  more  than  30  dB. 

Thus,  unless  something  is  done  to  limit  this  mechanism,  it  would  be  nearly  impossible  to 
construct  a  very-high-frequency  acoustic  microscope.  The  signal  would  deplete  too  rapidly  at 
higher  frequencies  (though  at  extremely  high  frequencies  the  dispersion  is  downward  at  all 
pressures,  as  shown  by  the  millimeter  path  lengths  of  10-20  A  phonons).33  Fortunately  in 
superfluid  helium  the  dispersion  is  adjustable.  At  low  pressures,  the  dispersion  is  upward,  and  so 
the  three-phonon  collision  occurs  at  a  finite  angle.  As  pressure  is  applied  the  dispersion  becomes 
more  linear  (collision  angle  becomes  smaller),  and  finally,  above  18  or  20  bar,  the  dispersion  turns 
downward.  With  the  dominant  scattering  process  thus  forbidden,  the  fluid  should  be  much  more 
transparent  for  sound  waves  at  helium  bath  pressures  above  20  bar,  and  the  microscope  should 
have  a  higher  signal-to-noise  ratio. 
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At  this  point  it  is  important  to  understand  the  difference  between  nonlinearity  and  dispersion. 


Nonlinearity  is  dependent  on  the  INTENSITY  of  the  beam  and  gives  the  relative  coupling  or 
interaction  between  different  modes.  Dispersion  is  dependent  on  the  frequency  of  the  mode,  and 
describes  whether  or  not  that  interaction  grows  as  the  beam  propagates  (or,  in  other  words, 
whether  momentum  is  conserved). 

2.  Experimental  Results  Using  a  Focused  Beam 

To  study  these  nonlinear  effects  at  high  power  densities  an  acoustic  lens  and  a  reflector  were 
placed  in  a  pressure  cell  (Fig.  3.2).  The  cell  has  been  described  previously.34  The  cell  allows  the 
lens-sample  spacing  to  be  adjusted  while  varying  the  pressure  applied  to  the  fluid  from  saturated 
vapor  pressure  (a  very  good  vacuum)  up  to  26  bar,  above  which  the  liquid  helium  freezes.  The 
experiment  measured  the  intensity  of  sound  pulses  which  bounce  off  the  reflector  (polished,  flat 
sapphire)  and  return  to  the  acoustic  lens.  The  received  signal  is  plotted  in  Fig.  3.3  as  a  function  of 
pressure  in  the  cell  for  two  frequencies  with  two  different  acoustic  lenses.  It  can  be  seen  that 
although  the  shape  of  the  curve  is  considerably  different  between  the  two  lenses  (believed  to  be  due 
to  slight  variations  in  the  shape  of  the  foci  caused  by  slight  variations  in  the  illumination  of  the 
lenses  by  the  transducers)  there  is  a  20  dB  increase  in  the  signal-to-noise  ratio  near  20  bar. 

A  study  was  also  made  of  the  effects  of  very  high  power  densities  at  various  pressures  using 
the  8  GHz  lens  and  flat.35  Plots  of  the  output  power  as  a  function  of  input  power  at  various 
pressures  are  shown  in  Fig.  3.4,  where  the  input  power  has  been  rescaled  in  terms  of  the  Mach 
number,  given  as: 


(3.1) 


where  po  and  Co  are  the  density  and  speed  of  sound,  respectively,  of  the  helium,  and  I  is  the 
intensity  of  the  sound  beam  (all  the  signal  power  is  assumed  to  fall  on  a  spot  which  is  diffraction- 
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FIGURE  3 . 2  Cutaway  view  of  the  pressure  cell,  showing 

(A)  Sapphire  rods  for  lens  and  flat  reflector,  (B)  and 
(Q  Miniature  bellows,  (D)  Pushrod,  (E)  Control 
capillaries,  and  (F)  Impedance-matching  circuitry, 
(from  Reference  12). 
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FIGURE  3 . 3  Improvement  of  the  Signal-to-Noise  ratio  upon 

pressurization  of  helium  at  frequencies  of  7  and  8  GHz. 
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FIGURE  3 . 4  Power  received  (0  dB  corresponds  to  the  noise  level  of 
system)  as  a  function  of  input  intensity,  given  as  Mach 
number  at  focus,  for  various  bath  pressures. 
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limited  with  diameter  of  one  wavelength  times  the  f-number  of  the  lens).  Three  regions  of  pressure 
with  different  behaviors  are  observed.  Each  of  these  three  regions  will  be  discussed  in  turn. 

At  low  pressures,  8  bar  and  12  bar,  the  dispersion  curve  is  upward,  or  "anomalous,"  that 
is,  the  speed  of  sound  increases  with  increasing  frequency.36  Spontaneous  decay  of  the  signal  is 
thus  possible  using  a  three-phonon  process  (3pp)  where  the  decay  products  have  frequencies  C0i, 
g>2  and  wavenumbers  kj,  k2  satisfying 

coi  +  0)2  =  cop;  kj  +  k2  =  kp  (3  2) 

where  the  subscript  p  (=pump)  refers  to  the  input  signal.  The  beam  can  thus  spontaneously  decay 
into  two  noncollinear  phonons.  Pump  depletion  is  the  result  of  parametric  amplification  of  the 
stimulated  decay  products.  Such  behavior  has  already  been  reported  for  plane  waves  in  superfluid 
helium.27 

At  an  intermediate  pressure,  15  bar,  the  dispersion  is  almost  linear.  As  the  dispersion  gets 
more  and  more  linear,  the  available  angle  between  the  decay  products  gets  smaller  and  smaller. 
Hence  the  sum  processes  should  dominate.  That  is,  the  phonons  in  the  beam  are  more  likely  to 
scatter  off  each  other  to  produce  harmonics  than  to  scatter  off  the  noise  to  decay.  The  beam  power 
thus  saturates.  The  power  remains  constant  over  a  15  dB  range  in  input  power. 

At  high  pressures,  19  bar,  20  bar,  22  bar,  the  dispersion  curve  is  downward,  so  that 
neither  the  three-phonon  decay  nor  harmonic  formation  is  phase-matched,  that  is,  conserving  both 
energy  and  momentum.  As  one  adds  increasing  amounts  of  power  to  the  system,  however,  the 
output,  rather  than  saturating,  continues  to  rise  to  a  much  higher  signal  level.  Then  as  more  power 
is  added,  the  output  level  actually  decreases.  The  pump  is  then  experiencing  more  than  the  simple 
saturation  that  might  be  expected  in  shock  wave  formation,  where  additional  energy  is  simply 
transferred  into  harmonics.  A  decrease  in  output  with  additional  input  would  seem  to  be  due  to  one 
of  two  causes.  Either  there  is  a  major  alteration  in  the  nature  of  the  medium,  or  there  is  a 
stimulated  process  occurring,  in  which  phase-matched  decay  products  are  parametrically  amplified 
at  the  expense  of  the  pump.  Each  will  be  discussed  in  turn. 


51 


3.  Analysis  of  Data 

A  number  of  other  effects  associated  with  high-intensity  sound  propagation  in  liquid  helium 
might  alter  the  medium  and  account  for  the  results  at  high  pressure,  namely  cavitation,  streaming, 
vorticity,  and  rapid  melting  and  freezing.  Above  a  certain  critical  amplitude,  cavitation  can  result  in 
massive  signal  loss.  There  was  not,  however,  any  increase  in  noise  associated  with  cavitation  and 
bubble  collapse,37’38  even  using  a  spectrum  analyzer.  Also,  evidence  from  another  system 
indicates  that  at  higher  frequencies  the  cavitation  amplitude  increases  linearly  with  frequency,38  and 
extrapolation  from  helium  experiments  imply  a  much  higher  threshold  than  even  the  intensities 
used  here  39  On  the  other  hand,  the  present  amplitudes  are  larger  than  the  onset  for  streaming 
reported  in  helium.40  Streaming,  however,  is  an  effect  usually  associated  with  continuous  sound, 
while  here  the  pulses  were  short  and  the  duty  cycles  low.  The  returned  signal  strengths  were  also 
independent  of  the  pulse  width  and  repetition  rates.  Furthermore,  observations  using  a  spectral 
analyzer  showed  no  spectral  (Doppler)  shifts  within  the  roughly  100  MHz  bandwidth  of  the  lens- 
amplifier  system.  Another  possible  effect  is  that  the  intensity  at  the  focus  could  be  high  enough  to 
cause  a  small  section  of  the  superfluid  to  go  "normal,"  due  to  the  high  vortex  density.  This 
process  should  be  a  slow  one.41  Experiments  using  multiple  pulses  showed,  however,  that  the 
pulses  were  completely  independent-noninterfering  on  time  scales  as  short  as  100  nanoseconds. 

A  final  alternate  explanation  is  rapid  melting  and  freezing  of  the  helium  at  the  focus.  This  process 
is  very  fast,  occurring  at  even  gigahertz  rates.42  However,  the  output  was  stable  in  time,  and 
showed  no  abrupt  change  with  increased  helium  bath  pressure,  as  might  be  expected  from  a  phase 
change.  Rather,  it  showed  a  slow  progression  between  the  various  regimes  of  pressure. 

Since  there  are  reasons  to  discount  the  helium-altering  effects,  the  possibility  of  a  stimulated 
process  must  be  considered.  At  high  pressures  the  simple  three-phonon  process  above  cannot 
conserve  both  energy  and  momentum.  The  next  simplest  decay  process  which  can  conserve  both 
is  the  four  phonon  process  in  which  two  pump  phonons  scatter 

coi  +  0)2  =  2c0p;  ki  +  k2  =  2kp  (3.3) 
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Referring  back  to  Fig.  3.4,  one  also  notices  that  the  rate  of  drop  of  output  signal  with  increasing 
input  signal  is  twice  as  steep  as  that  for  low  pressures.  This  result  is  consistent  with  Eqs.  (3.2) 
and  (3.3),  since  now  two  pump  phonons  are  consumed  in  each  stimulated  decay  event,  versus  one 
for  the  three-phonon  process.  Spontaneous  decay  (also  called  parametric  noise )  via  four-photon 
processes  has  been  observed  in  optics.43-44  We  believe  these  data  to  show  the  first  observation  of 
such  a  phenomenon  in  acoustics. 

This  conclusion  regarding  the  interpretation  of  the  collinear  gain  observed  in  the  parametric 
amplifier  should  be  adjusted  slightly.  The  probable  interaction  is,  like  that  one,  between  four 
phonons,  but  as  a  two-step  process  of  two  non-momentum-conserving  three-phonon  processes. 

IV.  The  Parametric  Amplifier  Experiment 

Studies  of  nonlinear  acoustic  interactions  in  superfluid  helium  at  temperatures  below  0.2*  K 
have  culminated  in  the  construction  of  an  all-acoustic  parametric  amplifier  at  gigahertz  frequencies. 
This  device  represents  the  shortest  wavelength  parametric  amplifier  ever  made,  with  signal 
wavelengths  shorter  than  1000  Angstroms  and  pump  wavelengths  shorter  than  600  Angstroms.  In 
the  experiment,  plane  waves  are  mixed  at  a  predefined  angle  in  the  superfluid  helium.  Two  gain 
regimes  are  observed.  The  first  regime  is  a  noncollinear  phase-matched  process,  in  which  the 
signal  phonon  stimulates  decay  of  the  pump  phonon  to  create  gain  at  the  signal  frequency.  The 
second  regime  is  a  four-phonon  collinear  process,  in  which  gain  on  the  signal  is  created  by  the 
second-harmonic  of  the  pump.  This  four-phonon  process  is  unusual  and  is  shown  to  be  a 
combination  of  three-phonon  processes  wherein  the  lack  of  conservation  of  energy  and  momentum 
in  one  process  is  mostly  compensated  by  the  other. 

1.  The  Plan 

The  experimental  geometry  is  depicted  in  Fig.  4.1.  The  pump  and  signal  phonons  are 
launched  from  a  sapphire  substrate  into  liquid  helium  at  an  angle  determined  by  the  dispersion  in 
helium  and  the  frequencies  of  the  pump  and  signal  waves.  The  signal  phonons  are  amplified  and 
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FIGURE  4 .  1  Experimental  geometry  of  Parametric  Amplifier 

experiment  in  superfluid  helium.  The  "bottom  crystal" 
(right)  reflects  the  beam  back,  or  else  a  third  transducer 
can  be  deposited  on  its  backside. 


5A 


detected,  either  (after  reflection)  by  the  transmitting  transducer,  or  by  a  third  transducer  on  the 
other  flat 

To  determine  which  frequency  pairs  (pump  and  signal)  will  be  phase-matched,  a  calculation 
is  done  using  the  data  of  Rugar  and  Foster,  who  characterize  the  dispersion  in  terms  of  the 
parameter  a,  where: 

2 

£o(k)  =  ck  ( 1  +  ock  )  (4.1) 

with  0)  the  circular  frequency  (2ttf)  and  k  the  wavevector  (2tiA).  As  a  function  of  the  pressure  P 
(in  atmospheres)  their  data  form  a  straight  line  (r2  =  .99): 

a(P)  =  (1.531  -  0.1039  P(in  bar) )  *  10-20m2  f4  2) 


Referring  again  to  Fig.  4.1,  the  30°  bevel  will  cause  sound  to  refract  to  an  angle  (with  the  normal) 
of  roughly  0.6°  (actually  the  angle  will  also  change  slightly  with  pressure  due  to  the  change  in 
sound  speed  c).  Using  the  dispersion  relations  in  Eqs.  (4.1)  and  (4.2),  then,  for  example,  at 
2.6  bar,  a  4.5  GHz  pump  will  be  correctly  phase-matched  to  a  2.6  GHz  signal. 

2.  First  Experimental  Design 

The  first  experimental  attempt  at  constructing  an  all-acoustic  parametric  amplifier  in 
superfluid  helium  used  a  multifaceted-crystal  design  (see  Fig.  4.2).  The  transducers  are 
positioned  to  overlap  as  the  beams  exit  the  sapphire  piece  at  the  bottom,  near  the  bevel 
intersection.  A  small  part  of  the  signal  from  the  transducer  on  the  top  bevel  (tilted)  surface 
reflects  from  the  bottom  bevel  (carefully  polished  to  be  parallel  to  the  upper  bevel).  It  was 
intended  that  this  return  could  be  used  to  tune  the  bevel  transducer,  so  as  to  optimize  its 
power  output  near  the  pump  frequency,  in  the  hope  of  maximizing  the  gain  of  the  amplifier. 

The  experiment  was  attempted  several  times,  but  the  reflected  signal  from  the  bottom  flat  (the 
helium  signal)  was  never  seen.  An  attempt  was  made  to  place  another  transducer  at  the  bottom  of 
the  cell,  to  operate  the  experiment  in  transmission.  The  shape  of  the  top  crystal  made  it  extremely 
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TRANSDUCER  SIDE 


PUMP  SIGNAL 

(BEVEL)  (MAIN) 


HELIUM  SIDE 


FIGURE  A .  2  Parametric  Amplifier  crystal  geometry  for  first 

experimental  try,  using  doubly-beveled  crystal.  The  flat 
bottom  (reflecting)  crystal  is  not  shown. 
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difficult  to  attempt  to  align  laterally  the  transducers  on  the  two  crystals,  so  the  system  could  not  be 
tested  in  water  (at  room  temperature).  Also,  despite  obtaining  return  sapphire  echoes  from  the 
bevel,  one  could  never  be  sure  just  how  well  the  beams  overlapped,  if  at  all. 

In  retrospect,  the  main  reason  for  failure  was  the  geometry  of  the  top  crystal.  Crystals  used 
in  superfluid  helium  acoustics  usually  use  a  matching  (antireflection  coating)  layer,  typically  of 
amorphous  carbon  with  a  thickness  of  a  quarter-wavelength  (for  the  center  frequency  desired). 

The  bevel  surface,  however,  meant  that  the  sapphire  surface  was  not  perfectly  flat  and  parallel  all 
the  way  to  the  edge,  since  polishing  techniques,  despite  great  care,  cannot  prevent  cracking  and 
splintering  close  to  the  edge.  The  carbon  matching  layer  thus  also  went  down  poorly,  with 
imperfections  near  the  edge  (which  is  precisely  where  the  experimental  beam  was  centered). 
Although  the  crystal  used  in  the  experiments  looked  visually  smooth  near  the  beam  exit  point, 
problems  in  adhesion  on  some  of  the  other  crystals  made  in  the  same  run,  in  some  cases  even 
peeling  off  the  carbon  layer,  suggest  that  even  the  experimental  crystal  may  have  had  enough 
unevenness  in  the  surface  and/or  matching  layer  thickness  (especially  on  the  1000  A  scale  of  the 
helium  wavelength)  so  as  to  give  a  poor  transfer  into  the  helium,  allowing  at  best  a  wave  of  poorly- 
defined  phase-front  to  be  transmitted. 

3.  New  Design 

The  new  design  adopted  when  the  decision  was  made  to  repeat  the  parametric  amplifier 
experiment  solves  many  of  the  worst  problems  of  the  first  design.  The  new  geometry  has  the 
bevel  surface  on  the  transducer  side  only,  with  the  carbon  side  (helium  side)  flat.  By  giving  up  a 
chance  to  tune  the  bevel  transducer  much  has  been  gained.  The  absence  of  a  bottom  edge  means 
that  the  sapphire  can  be  polished  more  easily  and  more  flat,  and  the  amorphous  carbon  matching 
layer  can  be  deposited  more  evenly,  for  good  beam  throughput. 

The  new  design  also  required  a  change  in  the  photolithography  technique  used  to  define  the 
transducers  for  "liftoff'  on  the  crystal  surface  (in  "liftoff*  the  gold  is  deposited  everywhere  atop  the 
photoresist,  and  upon  use  of  a  solvent  for  the  photoresist,  peels  away  everywhere  except  where  the 
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photoresist  has  been  previously  removed).  Under  the  old  design  the  main  (flat)  transducer  was 
defined  in  a  projection  microscope,  its  distance  to  the  edge  measured  (using  a  reticuled  eyepiece  in 
the  microscope),  and  the  position  of  the  bevel  transducer  calculated  by  a  program  taking  into 
account  the  walk-off  in  the  (anisotropic)  sapphire.  This  new  position  was  then  obtained  by 
looking  through  the  microscope  and  translating  the  sample,  using  the  reticule  and  small  specks  of 
dust  as  milestones.  The  distance  traveled  was  usually  two  to  four  fields-of-view  of  the 
microscope,  which  presented  a  dangerously  easy  source  of  error.  In  the  new  design,  without  the 
bottom  bevel  to  test  the  reflection,  the  placement  of  the  bevel  transducer  had  to  be  correct  the  first 
time— and  especially  since  the  first  run  had  had  no  results  it  was  determined  that  the  new  run  should 
eliminate  as  many  sources  of  error  as  possible.  (Later  it  was  also  learned  that  the  position  of  the 
beveled  transducer  "footprint"  could  be  determined  directly  using  small  latex  spheres,  discussed 
later). 

The  new  photolithography  technique  uses  a  single  exposure  through  an  electron- beam- 
defined  contact  mask  (see  Fig.  4.3).  The  calculation  of  relative  placement  was  thus  done  one  time 
(though  rechecked  many  times,  of  course,  before  programming  the  e-beam  machine  to  write  the 
mask)  and  the  mask  for  both  transducers  located  relative  to  the  bevel  edge  using  alignment 
"windows"  defined  in  the  mask.  Naturally,  since  the  crystal  has  a  bevel,  the  mask  will  be  a  contact 
mask  only  on  the  flat  surface  and  a  projection  mask  on  the  tilted  surface.  For  that  reason  both  of 
the  holes  were  made  below  size  (and  in  fact  the  "bevel"  hole  was  made  elliptical  so  its  projection 
would  be  round).  Thus  after  alignment  and  exposure  in  the  Kasper  mask  aligner  the  photoresist 
was  partially  developed  to  reveal  a  small  mark.  The  crystal  is  then  put  under  the  projection 
microscope  and  a  sharp  round  hole  of  the  correct  size  is  exposed  on  top  of  the  mark.  The  crystal  is 
then  placed  in  a  special  tilted  holder  and  a  round  mark  exposed  on  the  bevel  also.  In  that  way  the 
holes  in  the  photoresist  are  correctly  sized  and  sharp,  with  steep  walls  to  make  the  lift-off  process 
work  better. 
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collimated  light 
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FIGURE  4 . 3  Photolithography  of  transducers  for  Parametric  Amplifier 
experiment,  (a)  Top  view  of  mask  showing  alignment 
windows  and  exposure  circles,  adjusted  to  counteract 
acoustic  beam  walkoff  in  crystal,  (b)  Side  view 
showing  initial  exposure. 
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4.  Experiment  One 

The  first  experiment  produced  some  startling  results.  Gain  at  the  correct  "phase-matching" 
was  very  small,  and  a  larger  gain  was  seen  in  a  collinear  interaction  at  a  different  frequency. 

At  the  start  of  the  experiment,  gain  was  noticed  when  the  signal  and  idler  frequencies  were 
close  to  each  other,  surprising  since  for  the  planned  noncollinear  interaction  such  frequency  pairs 
are  very  phase-mismatched.  Typical  results  are  shown  in  Fig.  4.4,  which  monitors  the  output 
power  at  the  3360  MHz  signal  as  the  3800  MHz  pump  power  is  increased.  Notice,  from  (b),  that 
the  gain  is  still  present  on  the  beta  pulse;  in  fact,  it  is  larger!  Thus  the  dominant  interaction  here 
must  be  between  waves  traveling  in  the  same  direction.  Figure  4.5  shows  the  maximum  gain  as 
the  pump  frequency  is  varied,  thus  giving  the  bandwidth  of  the  process. 

The  behavior  of  the  pump  signal  was  also  monitored  as  the  signal  power  was  varied.  Figure 
4.6  (a)  shows  the  power-out/power-in  curve  for  the  pump  alone  (no  signal  power).  The  graph 
shows  the  usual  linear  climb,  then  sudden  saturation  and  drop  as  the  result  of  runaway 
amplification  of  the  noise.  The  second  small  additional  increase  some  30  dB  after  the  saturation 
point  possibly  represents  the  detection  of  an  outer  lobe  of  the  bevel  signal  (since  the  bevel 
transducer  is  significantly  less  efficient).  Also  shown  are  two  points  marking  the  "small"  and 
"large"  powers  chosen  for  part  (b).  Figure  4.6  (b)  shows  the  behavior  of  a  fixed  input  pump 
power  as  the  signal  power  is  increased.  In  both  cases  the  pump  is  depleted,  but  note  that  the  small 
pump  actually  experiences  gain!  For  power  to  flow  to  the  higher  frequency,  it  must  be  the  case 
that  the  driving  "effective"  pump  is  the  harmonic  of  the  stronger  input  beam. 

After  some  search,  gain  was  finally  found  where  it  was  supposed  to  be.  Figure  4.7  shows 
the  maximum  gain  on  the  alpha  signal  (there  was  no  gain  on  the  beta  signal)  for  the  calculated 
noncollinear  phasematched  frequencies  of  2.545  GHz  for  the  signal  and  various  pump 
frequencies,  at  an  intermediate  pathlength,  d=0.72  mm.  More  details  on  all  of  this  can  be  found 
in  the  thesis  by  Wright.45 
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Power  Out  at  3360  MHz  Signal  (dB)  Power  Out  at  3360  MHz  Signal  (dB) 


(a)  Alpha  Pulse  Gain 


(b)  Beta  Pulse  Gain 


FIGURE  4 . 4  Signal  Power  Received  as  Pump  Power  is  varied  for  (a) 
alpha  and  (b)  beta  pulses.  Both  experience  gain. 
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FIGURE  4 .5  Maximum  Gain  of  alpha  pulse  of  signal  as  a  function  of 
pump  frequency,  thus  giving  Gain  Bandwidth. 
d=0.70  mm. 
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Output  Pump  Power,  alpha  pulse  (dB)  Output  Pump  Power,  alpha  pulse  (dB) 


(a)  3800  MHz  Pump  Only 


(b)  3360  MHz  Signal  Applied 


large  pump 
small  pump 


FIGURE  4 . 6  Behavior  of  Pump  in  Collinear  Gain  Process  (alpha 

pulse):  (a)  Power  Out  vs.  Power  In  for  Pump  alone  (no 
signal  power  applied),  (b)  Pump  Out  as  Signal  Power 
increased  for  two  different  input  pump  powers.  Note 
gain  on  small  pump. 
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Gain  at  2545  MHz  Signal  (dB) 


FIGURE  4 . 7  Noncollinear  gain  bandwidth,  d=0.72  mm.  Gain  is  on 
alpha  pulse  only  (no  gain  on  beta  pulse). 
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V.  Acoustic  Microscopy  at  15.3  GHz 

At  the  conclusion  of  this  program  we  furnished  a  multilayer  transducer  to  M.S.  Muha  at  the 
Aerospace  Corporation  for  their  use  in  a  Scanning  Acoustic  Microscope  at  15.3  GHz.  The 
following  is  abstracted  from  their  publication.46 

"The  highest-resolution  acoustic  microscopes  have  been  diffraction  limited,  so  to  achieve 
higher  resolution  an  acoustic  microscope  must  utilize  a  shorter  wavelength.  The  shortest 
wavelength  that  can  be  used  depends  on  the  choice  of  the  fluid  that  couples  the  acoustic 
energy  from  the  lens  to  the  object.  An  acoustic  microscope  using  water  as  the  coupling 
medium  achieved  a  resolution  of  2000  A  at  a  frequency  of  5  GHz.  Further  improvement 
in  resolution  for  a  water  acoustic  microscope  is  hampered  by  a  low  signal-to-noise  ratio, 
which  is  caused  by  an  increasing  acoustic  attenuation  in  water  as  the  frequency  is 
increased.  Better  resolution  can  be  achieved  by  using  as  the  coupling  fluid  unpressurized 
superfluid  4He  at  temperatures  below  0.2  K,  because  of  its  low  acoustic  attenuation  and 
low  sound  velocity.  An  acoustic  microscope  operating  in  unpressurized  superfluid  4He  at  a 
temperature  of  0.1  K  achieved  a  resolution  of  200  A  at  a  frequency  of  8  GHz.  Further 
improvement  in  resolution  in  such  a  cryogenic  acoustic  microscope  is  limited  by  a  low 
signal-to-noise  ratio,  which  results  when  acoustic  nonlinearities  in  the  superfluid  helium 
cause  the  received  signal  to  saturate  and  then  decrease  as  the  input  power  is  increased.  The 
maximum  received  signal  (at  saturation)  decreases  as  the  frequency  of  operation  is 
increased.  The  acoustic  nonlinearity  believed  to  be  responsible  for  this  limitation  in 
received  signal  is  the  three-phonon  spontaneous  decay  of  the  input  signal  at  a  frequency  of 
CDp  and  a  wavenumber  of  kp  into  decay  products  at  frequencies  coi  ,  0)2  and 
wavenumbers  ki,  k.2  ,  where  energy  and  momentum  conservation  must  be  satisfied,  i.e., 

0)i  +  C02  =  top ;  ki  +  k2  =  kp  (1) 
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This  spontaneous  decay  process  is  self-enhanced  in  a  parametric  process  at  high  sound 
intensities  such  as  are  found  at  the  focus  of  an  acoustic  lens. 

The  acoustic  microscope  reported  here  eliminates  parametrically  self-enhanced 
spontaneous  decay  by  employing  4He  pressurized  to  greater  than  20  bar.  For  Eq.  (1)  to  be 
satisfied,  it  is  necessaiy  for  the  sound-dispersion  curve  to  be  anomalous;  that  is,  the  speed 
of  sound  must  increase  with  increasing  frequency.  Below  20  bar,  4He  has  an  anomalous 
sound-dispersion  curve.  Pressurizing  the  superfluid  helium  to  above  20  bar  causes  the 
dispersion  curve  to  be  normal  (i.e.,  the  speed  of  sound  decreases  with  increasing 
frequency)  at  all  frequencies,  so  that  three-phonon  spontaneous  decay  [Eq.  (1)]  is 
forbidden.  In  fact,  the  spontaneous  decay  of  a  single  phonon  into  any  number  of  phonons 
is  prohibited.  The  received  signal  is  then  no  longer  limited  by  spontaneous  decay.  A 
different  saturation  mechanism  still  occurs  in  helium  pressurized  above  20  bar,  but  onlyu  at 
a  much  higher  power  level.  For  example,  D.  R.  Wright  et  al.,  measured  an  increase  of  15 
dB  in  maximum  received  signal  of  when  the  pressure  was  increased  from  8  to  22  bar. 

They  attributed  the  saturation  in  received  signal  in  helium  pressurized  above  20  bar  to  a 
four-phonon  process,  where  two  signal  phonons  combine  to  produce  two  new  phonons  at 
different  frequencies. 

Pressurizing  the  helium  above  20  bar  also  eliminates  another  source  of  acoustic 
attenuation:  the  scattering  of  the  signal  phonons  by  thermal  phonons  in  a  three-phonon 
process.  This  process  is  forbidden  by  conservation  of  energy  and  momentum  in  a  fluid 
having  normal  dispersion.  With  this  thermal -phonon  scattering  eliminated,  the  only 
significant  source  of  attenuation  in  the  fluid  below  0.9  K  is  the  four-phonon  scattering 
referred  to  above.  At  temperatures  above  0.9  K  the  attenuation  becomes  large,  possibly 
because  of  roton  scattering.  Hence  a  microscope  using  superfluid  4He  pressurized  above 
20  bar  can  be  operated  up  to  approximately  0.9  K,  while  an  unpressurized  microscope 
must  be  operated  below  0.2  K  to  avoid  thermal-phonon  scattering.  The  acoustic 
microscope  described  here  is  cooled  by  a  sorption-pumped  ^He  cryostat  to  0.4  K,  while 
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the  unpressurized  4He  microscope  reported  earlier  required  a  dilution  refrigerator  for 
cooling  below  0.2  K. 

The  heart  of  the  acoustic  microscope  is  a  c-axis-oriented  sapphire  rod  having  an 
acoustic  transducer  on  one  end  and  an  acoustic  lens  on  the  other.  The  acoustic  lens  element 
is  a  polished  spherical  depression  ground  on  one  end  of  the  sapphire  rod.  An  800-A-thick 
layer  of  amorphous  carbon  evaporated  on  the  lens  surface  serves  as  a  matching  layer  to 
reduce  the  large  acoustic  impedance  mismatch  between  sapphire  and  liquid  helium.  A  400 
A  thick  layer  of  molybdenum  is  evaporated  before  the  carbon  is  deposited  to  increase  the 
adhesion  between  the  carbon  film  and  the  sapphire  substrate  and  to  improve  the  impedance 
matching.  The  acoustic  transducer  is  a  multilayer  ZnO  piezoelectric  thin  film.  Each 
layer  is  1700  A  thick,  and  the  crystal  orientation  of  alternate  layers  is  varied  to  attain 
periodic  coefficients  of  piezoelectric  coupling.  This  type  of  multilayer  transducer  presents 
a  higher  impedance  to  the  incoming  microwave  pulses,  resulting  in  a  better  match  and  a 
higher  conversion  efficiency  than  would  be  possible  at  15.3  GHz  with  the  single-layer 
transducers  previously  used.  The  measured  one-way  conversion  loss  of  the  multilayer 
transducer  is  15  dB." 

This  instrument  is  now  operational.  It  is  used  for  recording  images  of  IC  devices  with  a 
resolution  of  150  A. 
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